Europhysics Letters 



PREPRINT 



Full counting statistics of incoherent Andreev transport 

W. Belzig 1 and P. Samuelsson 2 

1 Departement Physik und Astronomie, Universitdt Basel, Klingelbergstrasse 82, CH- 
4-056 Basel, Switzerland 

2 Departement de Physique Theorique, Universite de Geneve, CH-1211 Geneve 4, Switzer- 
land. 



PACS. 74.50.+r - Proximity effects, weak links, tunneling phenomena, and Josephson effects. 

PACS. 72.70.+m - Noise processes and phenomena. 

PACS. 73.23.-b - Electronic transport in mesoscopic systems. 



Abstract. - We study the full counting statistics of heterostructures consisting of normal 
metal parts connected to a superconducting terminal. Assuming that coherent superconduct- 
ing correlations are suppressed in the normal metals we show, using Keldysh-Nambu Green's 
function, that the system can be mapped onto a purely normal system with twice the number 
of elements. For a superconducting beam splitter with several normal terminals we obtain 
general results for the counting statistics. 



Introduction. - A complete statistical description of a transport process can be obtained 
from the full counting statistics (FCS) [Q. The FCS of charge transport in mesoscopic con- 
ductors has recently attracted a lot of attention [2] . It was shown that scattering between un- 
corrected Fermi leads is described by binomial statistics , in which charge is transfered in 
elementary units. Subsequently the FCS in normal-superconducting systems have been stud- 
ied and it was shown that Andreev reflection leads to a doubled charge transfer [IflSKI . Later 
on, several approaches to FCS have been put forward. A theory of full counting statistics 
based on Keldysh-Green's function was developed |7J. This formulation allows a straight- 
forward generalization to superconductors HHH, mu hi-terminal structures ^U] an d Coulomb 
blockade systems ^J. As another example, the FCS of charge transport was expressed via 
the counting statistics J2| of photons emitted from the conductor. Very recently, a classical 
path- integral approach to FCS was developed jT3| . 

The second moment of the FCS is related to the shot noise Suppression of the shot 
noise from the classical Poisson value can result from the Fermi statistics of the electrons |15l 
llt>U17| . Particularly interesting are the correlations between currents at different terminals - 
the cross correlations. For uncorrelated Fermi leads they are in general negative [18117) . which 
has been experimentally confirmed 19 20121) . If the current is injected from a superconducting 
terminal, the sign of the cross-correlations may become positive [22 23, 24 25, 26 27]. If the 
proximity effect is absent both negative cross correlations |28I6| and positive cross correlations 
|2H] were theoretically predicted. 
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The positive correlations are a direct consequences of the charge being transferee! in pairs 
across the normal-superconducting interface 4 , 26, 2.3001 • This is the case, independent of 
the presence or absence of the proximity effect in the normal metal. Since the proximity 
effect introduces additional pair correlations, it is of special interest to investigate, under 
general circumstances, the charge transfer in normal-superconducting systems with suppressed 
proximity effect. 

Circuit Theory of Incoherent Andreev Transport. All transport properties of quasi- 
classical superconductor- normal metal heterostructures can be calculated from the so-called 
circuit theory p3U[32J|71 l 1 UU33 j . This theory is a discretization of the normal metal parts into 
connectors and nodes. The mesh has to be chosen in a way that approximates the real struc- 
ture to the required precision. To each node (label i) a matrix Green's function Gi is ascribed. 
A matrix current 7y flows through the connector between nodes i and j . Conservation of the 
matrix current on each node 7y = together with the normalization condition Gf = 1 
determine the unique solution for the whole circuit. Connections to the external world are 
represented by terminals (labeled with a) with a fixed matrix Green's function G Q . 

Electrons and holes (connected by Andreev reflection at the superconducting terminals) 
propagate through the structure. In the circuit theory the decoherence due to the finite energy 
difference 2E between electron and hole is described by a decoherence terminal d, which is 
connected to the respective node |32| . The matrix current from the node into the decoherence 
terminal has the form = — i(e 2 /h)(E/S)[a z , Gj\, where 5 is the spacing. Alternatively the 
coherence is suppressed by a magnetic flux <I> penetrating the node. The corresponding matrix 
current reads Ii ~ e 2 ($/$o) 2 fezGjO^, Gi] . Here a z denotes the third Pauli matrix in the 
Nambu matrix space and $0 is the flux quantum. 

The full matrix current conservation on node i then takes the from 



^4 + / id = o. (i) 



If we assume that the relevant energies are large in comparison to the mean level spacing 5 
multiplied by the dimensionless conductances of the connectors to the node (or if the magnetic 
field is large), all off-diagonal elements of Gi in the Nambu space are suppressed. To show this, 
we choose a representation of the Green's functions, in which the electron- and hole-Keldysh 
Green's functions are block-diagonal. For example, a normal terminal (labeled a) has the 
form 

Ga ~{ Gl{ Xa ))- (2) 
Here electron (hole) Green's function are still matrices in Keldysh space, defined as 

(1 — 9 f e (' 1 ' — 9 f e W \ 

-2(1 - fl ih) ) 2f a {h) ~ 1 ) P) 

with fa (h) (E) = (1 + cxp[(£ =F eV^/ksT])- 1 . The Green's function of the node i is 



&-( % % )• 



The decoherence matrix current is 
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Fig. 1 - Mapping of a beam splitter circuit, a) A number of normal terminal (here: three) is 
connected to one node. The node is connected to a superconducting terminal by a contact with 
transmission eigenvalues {T n }. Strong decoherence of electrons and holes is accounted for by a 
decoherence terminal connected to the node. This allows to map the circuit on the structure depicted 
in b). All normal parts of the structure a) are doubled and the superconducting terminal connects 
these two circuits. The connector between electron and hole circuit has transmission eigenvalues {R^}, 
where = T%/(2 — T n ) 2 is the Andreev reflection probability corresponding to the transmission 
eigenvalues of the original circuit a). 



Thus, if the decoherence current dominates over all other currents, the matrix current conser- 
vation equation Eq. becomes lf ec = and the off-diagonal components of Gi vanish, i. e. 
Gi becomes block-diagonal. 

The connection between two nodes is described in general by a matrix current 

f 2e 2 ^ r„ [&,(?,-] 

/^4 + r n ({G,G 3 }-2)' W 

where {T n } is the ensemble of transmission eigenvalues. The matrix current between two 
nodes with block-diagonal Green's functions is also block diagonal. Each block has again the 
form Thus, in the whole circuit (except at the superconducting terminals) the electron 
and hole blocks decouple. 

At energies well below the superconducting gap, a superconducting terminal at zero chem- 
ical potential has the Green's function 

«-(?;)■ <7 » 

The matrix current to the superconducting terminal has the form (JBJ, where the Green's 
function on the normal side has a block-diagonal form. Then, the block-diagonal components 
of the current can be rewritten, after some algebra, as 

Here we introduced the Andreev reflection probabilities = T%/(2 — T n ) 2 . Thus, the su- 
perconductor matrix current constitutes a (2x2)-matrix current between electron and hole 
Green's functions of the node. Note, that the 'transmission probabilities' are given by the An- 
dreev reflection probabilities. We conclude, that the transmission properties between electron 
and hole blocks (mirrored at the superconductor) have exactly the form of a normal contact 
with the usual transmission probabilities T n replaced by R^. 

With this result, we are ready to present the mapping rules illustrated in Fig. 
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1. All normal elements, normal terminals and all connectors between normal parts of the 
circuit are doubled (electron and hole circuit). 

2. The superconducting terminals are the only connections between electron and hole cir- 
cuit. They play the role of normal connectors, in which the normal transmission eigen- 
values T n are replaced by the Andreev reflection probabilities = T%/(2 — T n ) 2 . 

3. A counting field = ±x is assigned to corresponding electron and hole terminals. 

4. The counting statistics is obtained from the solution of the normal circuit. 

An important general property of the extended circuit is electron-hole symmetry. In our 
case this is formally a consequence of the perfect symmetry of the extended circuit. Mathe- 
matically the symmetry relation reads 

a x Gi(E)a x = -GH-E) . (9) 

Below we will use this symmetry relation to obtain some general properties of beam splitters. 

We can draw a general conclusion about the parity of the transfered charge number through 
the superconductor. Due to the Andreev process this number should be even. Since we 
have replaced the superconducting circuit by a purely normal on, the even parity is not 
obvious anymore. However, for a two terminal device the FCS is 27r-periodic in the difference 
between the counting fields Xe — Xhi which, according to our rules, should be replaced by 2%. 
Thus, the CGF is now a 7r-periodic function of x- This argument can be generalized to an 
arbitrary number of terminals, if we are interested only in the charge transfer through the 
superconductor. The total transfered charge is obtained by taking all counting fields equal. 
Again the CGF depends only on differences between counting fields and, therefore, charge 
transfers between two electron (or hole) terminals are not counted. With the same argument 
as given above, it follows that the total CGF is invariant under a shift of all counting fields 
by 7T. 

Counting statistics of two-terminal devices. - As a check, we first apply the mapping 
rules to a simple two-terminal contact without internal structure. The results for the cumulant 
generating function (CGF) of a normal contact with transmission eigenvalues {T n } is 
(k B T = 0) 

^NN(xi ! X2)=M^ln[l + T„( e i te-x 2 )_ 1) ] i (10) 

n 

where we have introduced the number of attempts per channel M = 2eVto/h. Using the 
recipe outlined above we replace T n by and xi ~ Xi by 2% and find 

S S N(x)=M£lii[l + itf(e i2 *-l)] , (11) 

n 

which is the result first obtained in Ref. 0|. 

The counting statistics of a diffusive wire between a normal and a superconducting terminal 
can also be found analytically. Let us write the CGF of a normal diffusive conductor with 
conductance gw as |7l l3*5ll3^] 

SmOa.Xa) = ^^acosh (2e ( * - 1) . (12) 

Now, according to the mapping rules, we have to replace the diffusive connector by two 
diffusive wires in series. We can neglect the interface resistance and use that the CGF for a 
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series of two diffusive conductors is the same as for a single diffusive conductor, 
conductance is halved and we thus find 



Ssn(x) 



9NVt 



acosh (2e l2x - l) 



The total 



(13) 



This result proves that the FCS of the diffusive SN-wire in the incoherent regime is the same 
as in the coherent regime |33| . which has so far only been demonstrated for the first two 
moments pOEHHl- 

As another example we consider a chaotic dot, which is much stronger coupled to the 
superconductor than to the normal terminal. In this case, the connector between the electron- 
and hole node does not contribute to the FCS (the Green's functions in the electron and the 
hole dots are the same). We can then consider the equivalent circuit with a single dot, 
symmetrically coupled to the electron and hole terminals. The Green's functions of electron 
and hole terminal are (for \E\ < eV and kgT = 0) 



G e W = T a g -(& x ±ia y )e Hx ' 



(14) 



For a moment we take independent counting fields in the electron and hole terminal, and 
replace Xe(h) ~ * ix on ly m the end. The solution for the central node is 



where Z = ^/exp(ix e — iXhj- Using this result, we find the CGF 



S( X e,Xh) = MJ2 1 * 



T 

- 1 - n 



1 + -y (yexp(zxe - ixh) - 1 



At this stage we can safely take the limit Xe(h) = ^X an d obtain 



S( Xe , Xh ) = Mj2^ [l + y( 



= i(xmod-7r) 



(15) 



(16) 



(17) 



which ensures the 7r-periodicity of the CGF. 

We now turn to the general case of the counting statistics of a beam splitter as depicted 
in Fig. The central node is connected to several normal terminals (labeled with a) and one 
superconducting terminal. This structure has been studied in different limits in Refs. |25ll2(:il 
I2"71l2"§] , From Eq. JSJl it is clear that the Andreev conductance gA = (4e 2 //i) ^ n i?^ governs 
the coupling between the electron and hole circuit. 

Counting statistics of a beam splitter - weakly coupled superconductor. We assume 
here that the superconductor is only weakly coupled to the beam splitter in the sense that 
gA -C <7e = ^2 a g a , where g a = (2e 2 /h) J2 n T n is the conductance of the connector to terminal 
a. All normal terminals are held at the same potential V and zero temperature. Expanding 
all quantities to first order in ^a/ss, we find (similar to Ref. |2§]) that the total CGF can be 
expressed by the CGF of the superconducting contact alone: 



S = M^Trln 



1 



4 



({' 



GeO ±. /^iev ~ 
,a x Lf c a. 



eO 



(18) 



Here is the electron Green's function of the central node in the absence of the supercon- 
ductor (for gA = 0). All Green's functions are evaluated at E < V and we used Eq. J3J to 
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express S in terms of Gf only. The Green's function G^° can be obtained quite generally. 
Due to the triangular shape of all G e a (see Eq. {2}), also G e ® has the same form. The matrix 
current between the central node and terminal a then becomes I a — ^-[G^ ,G%]. Thus, all 
normal connectors behave as tunnel contacts. The solution for the central node is j2E] 

Gf = —Y^g a G e a = a z + {a x + ia v )K, (19) 

where A = ^2, a Va& %XoL with Pa — ffa/ffs- It then follows from Eq. (|18fl . that the CGF is 

S({ Xa })=Mj2^[l + R*(A 2 -l)} . (20) 

n 

Such a CGF leads to a FCS of the form 

P( Nl ,N 2 , ...) = Pc(2Q)\ (l + (-1) 2Q ) jJggji- P *p»> ■ ■ ■ (21) 

where Q — ^2 a N a /2 is the number of transfered Cooper pairs. The total probability dis- 
tribution is therefore the probability Pc(2Q) that Q Cooper pairs are transfered, and then 
distributed among the normal terminals with respective probabilities p a . A similar result has 
previously been obtained for tunnel contacts and a chaotic cavity • Here, this deriva- 
tion holds for any type of connector, i. e. point contacts, diffusive wires, single transparency, 
etc. 

We note that if the beam splitter is connected to the superconductor by another connector 
(incoherent!), the solution is again given by Eq. I|2U|I . This is the case since the additional 
node, connected directly to the superconductor, has the same Green's function as the central 
beamsplitter node, Gf , to leading order in i- e. with the superconductor completely 

decoupled. 

Counting statistics of a beam splitters - strongly coupled superconductor. Another 
important case is when the superconductor is strongly coupled to the beam splitter, i. e. 
gA^> In this limit the superconducting connector is absent and the system is a chaotic 
dot coupled to the electron and hole terminals. We consider here the case, in which the normal 
terminals are coupled by tunnel contacts. The solution for the central node is then |26| 

6 e = K/V¥ , K = J2 9 y(gI + g£) . (22) 



The CGF follows straightforwardly, 

S({ Xa }) - M (TrV^ - se ) = M 



g a e l *<> - 



(23) 



In the present form the 7r-periodicity is evident. However, calculating explicitly the cumulants 
by taking successive derivatives of S({x a }) with respect to Xa, we can equally well take the 
square-root in Eq. (J53J), giving S({x a }) = MJ2 a g a [exp(ix a ) - 1] = J2 a S a(Xa)- From this 
form of S, no longer evidently pi-periodic, it follows, interestingly, that all cross correlators 
vanish. 
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Conclusions. - We have studied the full counting statistics of normal metal-superconductor 
heterostructures in the incoherent regime. The original circuit with one superconducting ter- 
minal can be mapped on a purely normal circuit consisting of an electron and a hole block. The 
superconductor plays the role of a normal connector between the electron and the hole block, 
with the usual transmission probabilities {T n } replaced by the Andreev reflection probabilities 
{R^ = T%/(2 — T n ) 2 }. We have illustrated our approach with several examples. 

* * * 
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